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This is a status report on where PHENIX stands in terms of mapping out the landscape 
of jet correlation in pt, hadron species, i/s. We discuss separately high pt correlation 
results and low and intermediate px correlation results. The former is sensitive to the 
quenching of the jet by medium, the later allows detailed study of the medium response 
to the (di-)jct. 
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1. Introduction 

One of the primary focus of the RHIC program is to study the properties of the 
strongly interacting quark-gluon plasma (sQGP) with high transverse momentum 
(Pt) jet and dijets. The properties of the plasma can be deduced in two ways: on the 
one hand, the suppression of single jets and di-jets signals at high pt tells us about 
the stopping power or the transport properties of the medium. On the other hand, 
study the the response of the medium to the jets, especially their energy dissipation 
in the medium, can reveal the collective properties or equilibration process of the 
medium. 

The studies of two particle azimuth correlation at RHIC have revealed detailed 
information on the jet interaction with the medium, beyond what we learned from 
the single particle spectra measurements. Experimentally, any deviations of the 
correlation pattern in heavy-ion collisions from that in baseline p+p collisions are 
attributed to some sort of medium effects. Such deviations are found to be depen- 
dent on the pt of the trigger hadrons (pt,t) and partner hadrons (pt.a)- FigfT] 
show a schematic sketch of the landscape of the two particle correlations in pt,t 
or pt,a- Based on existing experimental results and theoretical calculations, the 
while pt range is divided into four regions, each of which is probably sensitive to 
different physics. At high pt, STAR HI shows that the away side jet magnitude is 
suppressed by factor of 4-5 relative to p+p, but maintaining a relatively unmodified 
shape. It is consistent with the fragmentation of the primary jets that survives the 
medium either through tangential emission ^ or punching-through mechanisms ^ 
or both S At intermediate pt range(l-4 GeV/c), the away side jet is broadened El 
or even concave like ^ in shape but its amplitude is enhanced relative to p+p 
The features were believe to be the consequence of the response of the medium to 
the energy degradation of the high pt jets. Many theoretical scenarios have been 
proposed in the last few years to describe the e nergy dissipation processes, includ- 
ing gluon radiation ^ Cherenkov radiation^HE "Mach Shock" mechanism El. But 
deflected jet pictures such as hydrodynamical wakeQ^or multiple scattering of the 
jet by the medium ^Sl W ere proposed as well. The first two models require addi- 
tional particle production, whereas the other models require redistribution of the 
momentum/energy among existing particles. 

In between these two pt regions, experimental data show an suppressed and 
possibly very flat away side with no visible peak structure This region (we call 
moderately high pt region) could have contributions from both jet fragmentation 
and responses of the medium. The flat but still positive structure could be the 
combined result of a fragmentation contribution concentrated around A(f> = -k and a 
medium contribution centered at 1 radian away from n. The flat structure probably 
implies the two components have roughly equal yield at this pt region. 

Lastly is the low pt region. In most jet quenching models USE the 
bremsstrahlung radiation was believed to be largely contained in a narrow an- 
gular range (0.5 rad) around the jet axis. These gluons might be the precursor of 
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Fig. 1. A sketch of possible contributions to the away side jet yield as function of pt of triggers 
and/or partners. The pt axis is schematically divided into four different regions. It is chosen such 
that the pt,t > PT,A (see Fig |16t ). 



the collective medium response, but it is also conceivable that some gluons simply 
thermalized with the medium after additional scattering. In this picture, one would 
also expect an enhancement in the region around n at low pt and it should have 
an exponential thermal shape. 

Given that the correlation patterns depend strongly on the pt of both particles 
in the pairs, it is crucial to map out the awayside jet properties in a broad pr range 
and to identify the transition regions between different correlation patterns. We 
focus separately on the high pt region, which is ideal for studying jet quenching 
and jet tomography, and on intermediate pt region, which is dominated by the 
medium response. PHENIX results on hadron-hadron, identified hadron-hadron 
and three particle correlation are discussed in this context, with an emphasize on 
new results we have shown in QM2006. 



2. px evolution of correlation pattern 

In correlation analysis, one correlate hadrons in one px region ( "trigger" ) with those 
in another pt window ("partner"). The hadron pairs from same jet tend to appear 
at A<p — \<pA — 4>b | ~ and those from for back-to-back dijet tend to appear at 
A(f> ~ 7r. We are interested in the partner yield distribution per-trigger, Yj et (PTY ), 
which is defined in a way similar to the previous d + Au and p + p analysiJ^H^ 



Y- Kt = + . V ;\ - £ (l + 2vZvS cos 2Ar 

J 2TrN t e \N ml *(Acf)) ^ V 22 S 

with the additional 1v\v^ cos2A0 term to take into account the flow modulation 
of the background pairs in Au+Au collisions. N t is the number of trigger, e is 
the efficiency for associated hadrons in full azimuth and in ±0.35 pseudo-rapidity. 
N is (A(j)) and 7V mlx (A0) are same-event pair and mixed-event pair distributions, 



February 8, 2008 21:56 WSPC/INSTRUCTION FILE jjia'iwcf06 



4 Jiangyong Jia 



respectively. The superscript t and a stands for the trigger and partner particles. £ 
is a normalization factor which is the ratio of the combinatorics pairs in the same 
event to those in the mixed event. £ is typically bigger but very close to 1. 

Fig[2]shows the PTY for various combination of trigger and partner px in central 
Au+Au collisions, in comparison with p+p. The partner px is constrained to be 
less than trigger px- From left to right and top to bottom, the px ranges of one or 
both particles increase, going from 3 — 4 x 0.4 — 1 GeV/c (trigger px vs partner px) 
to 5 — 10 x 5 — 10 GeV/c. The solid lines around the points represent the error due 
to the elliptic flow uncertainty, which peaks at and 7r and is the dominating error 
at px < 3 GeV/c. The shaded bands indicate the error due to uncertainty on the 
background level which is fixed by the Zero Yield At Minimum (ZYAM) method^^. 
This error depends on the statistical uncertainty at the ZYAM minimum and is the 
dominating at px > 3 GeV/c. 




Fig. 2. The yield of associated hadron per-trigger in A<j> for successively increasing px,T and Pt,A 
in 0-20 % Au+Au collisions. The two lines around the data indicate the errors due to uncertainty 
of elliptic flow; the shaded bands around indicate the uncertainty from ZYAM method. 



The top row shows the jet yield in three successive partner px ranges covering 
0.4 to 3 GeV/c with trigger px fixed in 3-4 GeV/c. The away side jet shapes 
are dramatic different between Au+Au and p+p. The p+p data always peak at 
7r; The Au+Au data show a concave shape with a dip around 7r and two side 
peaks around 1 radian from n. The concaveness grows with increasing partner 
px, however the location of the side peaks is not changing. To characterize the 
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distorted away side jet shape, we divided the away side into a "Head" region (tt/6 < 
| A.0 — 7r| ) and a "Shoulder" region (tt/6 < |A</> — tt\ < ir/2), as illustrated in Fig.lc. 
They are so named because for normal jet as in p+p, the "Head" region contains 
bulk of the jet fragmentation, whereas the "Shoulder" contains the tail of the jet 
fragmentation and the radiative contribution. The "Head" region is sensitive to 
the level of jet quenching while the "Shoulder" region is suitable for studying the 
medium response. Clearly, the behavior of the Au+Au data in the two regions 
are dramatically different from that of p+p. Namely there a suppression in the 
"Head" region and an enhancement in the "Shoulder" region. The suppression in 
the "Head" region sets in around 1 GeV/c and grows with larger partner px, while 
the enhancement in the "Shoulder" region persists up to 4 GeV/c. 

In the bottom panels of FigfSJ as the px of both hadrons are further increased, 
the awayside "displaced" peaks in Au+Au data seem to be compressed relative 
to p+p and the Au+Au near side. This reflects the fact that the away side yield 
drops faster than the near side with increasing partner px- In 4-5 GeV/c bin, the 
Au+Au away side is broader than p+p but is no concave. At the highest px bin 
(5-10 GeV/c), the away side turns into a convex shape, similar in shape to p+p but 
it's magnitude is largely suppressed. 

Note that medium modification is on the pairs. Thus the shape of the away side 
should be symmetric with respect to the trigger px and partner px- Figj2] suggests 
that the away side can be categorized in following four regions depends on the px,T 
and px,A based on its shape: 1) a flat region at px,T or px,A < 1 GeV/ c; 2) a concave 
shape region at 1 < px,t,Pt,a 4 GeV/c; 3) a convex region px,t,Pt,a ^ 5 GeV/c; 
4) a transition region (almost flat away side) Px,t{pt,a) ^ 5 GeV/c px,a(pt,t) ^ 3 
GeV/c. 

3. Medium response 

3.1. Hadron-hadron correlation 

PHENIX have tried different methods to quantify the behavior of the away side 
concave shape at intermediate px- The location of the side peaks are characterized 
by "D" Si the distance of the peak to 7r. "D" is found by triple gauss function fit 
of the away side: 

J(A0) = G N {A(f>) + G A (Acj) -tt-D) + G A (Acj) - tt + D) (1) 

Fig|3] show D for several collision systems and collisions energies. The position of 
the peaks was found to be around 1 radian from tt independent of the centrality 
( Apart >100), collision energy and system size. In addition, "D" was found to be 
relatively insensitive to the change of partner px ®, which is also confirmed in FigJ21 
This suggests that "D" reflects some kind of universal properties of the medium, 
such as the speed of sound as suggested by the M ach C one mechanism. It can't 
be described by Cherenkov gluon radiation models or most of the deflected 
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jet models such as hydrodynamical wake^lor Markov random scattering Such 
mechanisms would predict narrowing of the peak angle with increasing px- 



T3 
(0 



Q 0.5 



1 



J i p 



100 



★ d+Au 200 GeV 
■ Au+Au 200 GeV 
a Cu+Cu 200 GeV 

• Au+Au 62.4 GeV 
▼ Cu+Cu 62.4 GeV 

_L 



200 
N 



300 



400 



part 



Fig. 3. Collision centrality, energy, and system size dependence of shape parameters "D" 



PHENIX also try to characterize the convexity of the away side jet shape with 
the ratio of average jet yield in the "Head" region to that in the "Shoulder" region, 
Rm, 

_ / Head dA^ et (A0) / / Shouldcr dA^ ot (A^) 
HS Head / Shoulder [ ' 

Each of the distribution in Figj2] produces one value of i?ns that captures essence 
the the away side jet shape. In general, one expect i?us < 1 f° r a concave shape, 
Rns ~ 1 for a flat distribution and -Rhs > 1 for a convex shape. Thus Rhs is an 
ideal quantity for studying the transition between the "Shoulder" dominated region 
to "Head" dominated region in px- Note i?ns is purely a shape variable, thus it is 
symmetric with respect to px.x and px.A, i.e. Rhs(j>t,t,Pt,a) = Rhs(pt,a,Pt,t)- 

Fig[4] summarize the px,T and px,A dependence of i?ns i n P+P and central 
Au+Au collisions. It is presented in fine bins of partner pt for three trigger px 
regions. The wide red shaded error bars represent the elliptic flow error, which is 
correlated in px] The narrow green shaded bar represent the ZYAM error, which can 
vary from point to point. In p+p collisions, i?ns is always above one and increases 
with pt,a- This suggests that the away side is always peaked around it in p+p and 
its width narrows with increasing pt,a- The same ratio in 0-5% central Au+Au 
collisions are drastically different. For trigger pt in 2-3 GeV/c, i?ns decrease from 
around 1 at pt,a < 1 GeV/c to a level of 0.3 ± 0.1 at 3-4 GeV/c. For trigger pr 
in 4-5 GeV/c, the decrease of i?ns is less significant. The level at low partner px 
is around 1, relatively insensitive to the trigger px- This suggests that away side is 
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flat at low partner px independent of the trigger px- 
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Fig. 4. Rhs (The ratio of the average yield in head and shoulder region) as function of Pt,A in 
three Pt,t ranges for 200 GeV 0-5% central Au+Au and p+p. Concave i?HS < 1> na t ^?HS ~ 1> 
convex J?hs > 1- 

Fig(2]suggests that not only the away side jet shape but also the jet multiplicities 
are modified in central Au+Au collisions relative to that in p+p. The modifications 
depend on trigger px, partner px and A<f> and can be quantified by /aa, the ratio 
of PTY in Au+Au to that in p + p within a certain A0 window W: 

dA^ +Au / [ d^Yg* (3) 

The top panels Fig[5]show the Iaa as function of partner px at away side (| Ac/)— ir\ < 
it/2) for three trigger px bins. Jaa is greater than 1 at low partner px but drops 
towards high partner px- This implies that central Au+Au data have a significant 
enhancement at low partner px and a strong suppression at high partner px- For 
higher trigger px, the enhancement at low partner px is smaller and the suppression 
at high partner px is stronger. The point where Iaa cross 1 shifts to lower partner 
Px, indicating a stronger decrease of 7aa m partner p T for higher trigger p T - The 
observed px,x and px,A dependence is a consequence of the competition between the 
enhancement in the "Shoulder" region and the suppression in the "Head" region 
shown in Figj2] At low px,A, the away side is (thus Taa) is dominated by the 
"Shoulder" region; at sufficiently high px,A, the away side is dominated by the 
"Head" region. The point where Iaa = 1 would depend on both trigger and partner 
px- For 4-5 GeV/c trigger, the suppression level seems to saturates around 0.3 at 
Px,A > 3, close to the Raa value for 4-5 GeV/c inclusive hadrons. 

The bottom panels of Fig[2] show the Jaa as function of partner px at away 
side (\A(j)\ < 7r/3) for the same three trigger px bins. We see similar enhancement 
at low partner px, and the enhancement diminishes towards high partner px- As 
one increase the trigger px, the dependence on partner px weakens, and the Jaa 
distributions become flatter in partner px- The near side enhancement has been 



February 8, 2008 21:56 WSPC/INSTRUCTION FILE jjia'iwcf06 



8 Jiangyong Jia 



: 2.0<p <3.0GeV/c " 

1 1 r T,lrig 

1 I * 0-20% 

■ », . 


1 1 1 i 1 1 1 1 i 1 1 1 1 i 1 1 1 1 i 1 1 i 

: 3.0<p <4.0GeV/c " 

T.trlg 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

" 4.0<p <5.0 GeV/c " 

T.tng 


i 


.... I .... I .... I .... I . •. . 




) 1 2 3 4 E 12345 12345 

p (GeV/c) p (GeV/c) p (GeV/c) 

r T,assoc "T.assoc T,assoc 


2.0<p T <3.0GeV/c 

0-20°/o ,r ' 9 

I 1 '. 


-' ' ' ' ' ' ' ' 3.b<p T ' , <4.6Ge^/c' ' '- 

"T,trig 

Mi 

• 


' 4.b<p T <5.6 Ge^/c 


• — I — 

• 


■ 

• 

.... i .... i .... i .... i ... . 


* • • 

.... i .... i .... i .... i ... . 



n l .... i .... i .... i .... i .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I 

H) 1 2 3 4 E 12345 12345 



p (GeV/c) p (GeV/c) p (GeV/c) 

r T,assoc "T.assoc T,assoc 

Fig. 5. I\a for three different trigger bins in the away side (top panels) and the near side (bottom 
panels) in 0-20% central Au+Au collisions at 200 GeV. 



observed by STAR El and PHENIX Q2EI this enhancement was attributed by 
STAR to the a ridge contribution which extends out in Ary up to ±2. The ridge 
contribution seen by STAR in |A77|jl.7 is about factor of three compare to the 
pure jet contribution, and was argued to be reason for the near side enhancement. 
The enhancement seen in PHENIX is smaller than observed in STAR, because of 
the ridge yield in PHENIX's limited rj acceptance is smaller. The decrease of the 
enhancement with increasing trigger pt suggests that the ridge component is softer 
than the jet component. The px range where the ridge yield is important seems to be 
in Pt,TiPt,a < 4 GeV/c, very similar to the range for the enhancement at the away 
side. PHENIX have also observed a broadening of the near side jet width in central 
Au+Au collisions l2Ql at intermediate pr, the broadening was found to disappear at 
high trigger pt (also see Fig fT6|) . These observations are summarized in FiglBl 
The width results are also suggestive of a ridge contribution that is important 
at intermediate px but disappears at large px- Due to its limited r\ acceptance, 
PHENIX so far haven't be able to observe the ridge signature directly. However, 
the detailed study of near side jet shape and yield in broad pr range can provides 
indirect but still valuable constrains on the properties of the ridge. 

PHENIX have carried out the intermediate pr correlation analysis at both 
Y / snn = 200 and 62.4 GeV ®. CERES also carried out similar analysis at 
y/sNN = 17.2 GeV selection ^! The results from the three collision energies for 
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Fig. 6. Centrality dependence of the near side gauss width in 200 GeV Au+Au collisions at 
intermediate px of 2.5 — 4 X 2 — 3 GeV/c (Left) and at various high pj> bins above 5 GeV/c 
(Right). 



1 < Pt.a < 2.5 < pt,t < 4 GeV/c are summarized in the FigJT] The middle panel ^ 
shows the "jet fraction" , so the scale can't be compared directly with the other 
two panels which show the per-trigger yield. But the jet shape can be directly com- 
pared among the three distributions. FigjT] indicate a clear modification of away 
side jet at all energies. They all have displaced peaks around 1 radian from tt, the 
concaveness of the away side distributions seems to decrease for lower y/s. Since the 
trigger selection are similar, we expect the away side jet energy would be mostly 
determined by the trigger px, modulo a weak y/s dependence of the trigger bias 
(reflected by (z)) due to the y/s dependence of the parton spectra shape. CERES 
Time Projection Chamber rj range is 2.1-2.7 in lab frame, which corresponds to 
0.1-0.7 in CM frame. Thus its rapidity range is 0.6 and is close to PHENIX range 
of 0.7. Thus the jet yield at 17.3 GeV can be compared directly with that at 200 
GeV after multiplied by 0.7/0.6 = 1.17 to correct for the difference in r\ acceptance. 
The maximum and minimum value is 0.17 and 0.07 for 0-5% central Au+Au col- 
lisions at 200 GeV (FigEk) and 0.08 and 0.07 for 0-5% Pb+Pb collisions at 17.3 
GeV after the acceptance correlation (Fig[7t). Thus the amplitude of the displaced 
peak in CERES is about factor 2 lower than PHENIX value, whereas the level at 
the "Head" region in CERES is surprisingly close to PHENIX value. This proba- 
bly suggests that the yield in the "Head" region is dominated by fragments from 
jet. The jet multiplicity are somewhat higher at RHIC due to smaller (z), but jet 
quenching is stronger at RHIC than that in SPS. The yield at "Shoulder" region is 
significantly smaller at SPS than at RHIC, possibly suggesting a weaker medium 
effect. Clearly detailed study of the yfs dependence of the "Head" and "Shoulder" 
yield can provide crucial constrains on interplay between jet quenching and medium 
response. 

PHENIX also obtained some interesting but puzzling results for correlations 
between very low px hadrons (0.2 < px < 0.4 GeV/c) FigE] shows azimuthal 
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correlation in 0-5% Au+Au collisions for like-sign pairs (left) and unlike-sign pairs 
(right). The like-sign correlation shows a near side peak consistent with HBT cor- 
relation, whereas unlike-sign correlation does not. What is surprising is the ob- 
servation of a displaced peak at a location that is consistent with what has been 
observed for high px correlation. Whether it has the same physics origin is not 
clear. If this is due to correlation between soft particle in the jets or mini-jets, why 
the unlike-sign correlation does not show any enhancement at the near side as well 
(the ridge argument)? It is also worth pointing out that the displaced peak mostly 
concentrated in a very narrow window of |A?y| < 0.1 and that the peak only shows 
up in large A part > 100 and is insensitive to the collision energy. It is observed at 
0-5% Au+Au at 62.4 GeV but not in central Cu+Cu collisions at 200 or 62.4 GeV. 
Further characterization of the shape and yield of this away side peak is underway 
to clarify the picture. 
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Fig. 8. Low pt like-sign (left) and unlike-sign correlation function for 0-5% central 200 GeV 
Au+Au collisions. 
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3.2. Three particle correlation 

Three-particle correlation can provide direct topological features on the interaction 
of the jets with the medium, thus potentially allows clear distinction of various 
modification patterns such as Mach-cones and Deflected jets which can't be directly 
differentiated with two particle correlation method. There are four independent 
angular variables among the three particles, A^»i2 = 4>i — 02,A</>i3 = 4>i — (f>s, 
A?7i2 and A7713 . If one is mainly interested in the azimuthal correlation, as is done 
in STAR ^21 and in earlier PHENIX measurements ^Sl, then there are only two 
independent variable A4>i2 and A0i3. Recently PHENIX has chosen to characterize 
three particle correlations through the use a local polar coordinate system where 
the trigger direction is the z axis (see insert of Figj9]) ^£1 In this scheme, the four 
independent variables are B\fi\, <\>\ and 4>2- Note that the coordinate system is 
defined for each trigger so it changes trigger by trigger relative to the global frame. 
In the data analysis, PHENIX chose AO* = 0* — 6* 2 and Acf>* = </>* — $jj as the 
two independent variables. These two variables combine the information of all four 
angles, thus provide a very different way of presenting the jet signal. 




PHENIX Prdimiuaiy 



Fig. 9. (Left) Full A<f>* , A9* correlation surface in 10-20% Au+Au collisions. (Middle insert) 
Schematic illustration of coordinate system. (Right) the 1-D projection of along A0* for 1.65 < 
A9* < 2.2 rad. 

Left panel of Fig[9]shows the resulting three-particle correlation function in A9* 
and Acjf . The soft-soft and soft-hard flow backgrounds have been subtracted, but 
it still includes the 2+1 (two particles from di-jet and one from background) flow 
background. It is important to note that this correlation function as generated also 
contains the effect of the PHENIX acceptance. One can see three peaks in the ring 
region (1.65 < A8* < 2.2 rad) between the two gray circles, located at A<fi* w and 
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Acf>* « ±7r/2. These peaks were not observed for simulation where only normal di- 
jet is included, thus are signatures for the medium effects. After removing the 2+1 
flow contributions, right panel of Fig[9] shows the projection along the ring on the 
A<jf axis, together with expectation of simplistic cone jet and deflected jet scenarios 
given by Monte-Carlo simulations. In both scenarios, one expected a narrow peak 
around A0* « and a broad peak around Acfi* sa ir/2. The difference is that the 
deflected jet scenario predicts a bigger jet amplitude at A(f>* ss than that for the 
cone jet. From right panel of FigO it appears that the PHENIX measurements 
are most consistent with the excitation of a Mach cone in the medium and not 
consistent with a deflected jet hypothesis. 

3.3. Identified hadron-hadron correlation 

Previous results indicate that the jet properties at the nearside, awayside "Head" 
and awayside "Shoulder" regions are quite different. The near side is consist of 
a ridge component along Ai] and a jet fragmentation contribution. The awayside 
"Shoulder" region is dominated by collective medium response and the "Head" 
region seems to be consist of jet fragmentation plus possible feedback from the 
radiated gluons. Given the complexity of the problem and many theoretical possi- 
bilities, one need additional handles to test this picture. Correlation with identified 
particles can reveal the composition or the chemistry of the jet yield in the three 
regions, thus provide valuable constraints on the underlying physics. 

FigfTUl shows the partner meson and baryon yield in 1-1.3 GeV/c (top panels) 
and 1.6-2.0 GeV/c (bottom panels) for unidentified trigger in 2.5-4 GeV/c. Results 
for 0-20% and 20-40% Au+Au collisions are shown in the left panels and right pan- 
els, respectively. The partner baryon yields are scaled to approximately match the 
partner meson yields. Clearly, the away side shows a concave shape for both partner 
mesons and baryons. This is not surprising for partner meson since it dominates 
the charged hadrons. But it is interesting to see that partner baryons also have a 
similar but less concave shape (higher in the head region and lower in the shoulder 
region in the figure). 

Fig llll shows the ratios of the parter baryon/meson yield as function of partner 
Pt for several centrality bins at the near side (left panel) and away side (right panel) . 
The ratios are compared with that for the inclusive spectra and e + e~ collisions. 
The ratio for the near and away side jet grow steadily with increasing partner 
Pt in all centralities, similar to inclusive hadrons and normal jet fragmentation 
measured in e + e~ collisions. But the increase is a little bit stronger in central 
Au+Au collisions as well as for the away side. The p/(ir + k) ratio is about 0.1 for 
pure jet fragmentation and ~ 0.4 for inclusive spectra in central Au+Au collisions. 
The p/(tt + k) ratio for the away side jet is somewhere in between. If the Mach-cone 
scenario is true, it is conceivable that the jet excites the medium at the away side, 
which then fragments in the usual recombination prescription. This would naturally 
lead to the centrality dependent baryon/meson ratio in the away side. 
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Fig. 10. The associated meson and baryon A(f> distribution for unidentified charged hadron triggers. 
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Fig. 11. Ratio of partner baryon to meson at the near side (left panel) and away side (right panel). 
It is plotted as function of partner px for 0-20%, 20-40% and 70-90% centrality bin. Trigger 
particles are charged hadrons with 2.5 < Pt < 4.0 GeV/c. 



Fig[l2] shows the near side yield where the trigger is identified in 
PHENIX^left) compared with results from STAR (right) PHENIX Results 
show a steady increase with centrality, except a quick drop of the yield for trigger 
baryon at -/V par t >250. The increase trend is qualitatively consistent with results 
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Fig. 12. (left) the near PTY in \A<f>\ < 0.94rad for baryon-meson (squares) and meson-meson 
(circles) correlations as a function of N pa , r t- (right) Similar results from STAR. 



obtained from hadron-hadron correlation (see Fig|5]) , which was argued to be due to 
the ridge in Ary. Similar trend is also observed by STAR in the right panels, except 
that in their case, the yield for trigger baryon is significantly larger than for trigger 
mesons and yield for trigger baryon does not drop in central collisions. However we 
should note that the trigger px in STAR is different from PHENIX and they also 
include more ridge contribution due to a much larger Arj coverage. In any case, 
a good theoretical model should be able to simultaneously explain single particle 
yields, elliptic flow and as well as these PID correlation results. Quark recombina- 
tion (thermal-thermal and thermal-shower) was quite successful in describing the 
spectra and elliptic flow, but so far was not very su cces sful in describing the jet 
correlation yield, especially the centrality dependence L ^ JJ . 



3.4. Comments on V2 background subtraction 

One of the important systematic errors in the jet correlation comes from elliptic flow 
background subtraction (see EqQ]). In the past years, many different methods have 
been developed to measure the V2- Reaction-Plane (RP) Method (i>2{RP}), 2 par- 
ticle commulant method («2{2}), 4 particle commulant method (^{4}), and SMD- 
ZDC method (Reaction-plane for direct flow These measurement s have quite 
different sensitivities on the non-flow effects and eccentricity 

fluctuation™™. 

The obvious question to ask is what vi one should use in the jet correlation analysis? 

The two particle correlation method and two particle commulant method are 
closely related. Naturally, all non-flow effects and event by event fluctuations that 
affects the commulant v% would also contribute to the two particle azimuthal cor- 
relation. Thus it seems that «2{2} should be used in the background subtraction 
except that one would like to remove the jet bias to f2{2} since it is the signal. 
PHENIX use the W2{RP} in all correlation analysis. The U2{RP} is determined by 
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the BBC which sits in the forward region (3 < |?7| < 4). The RP method measures 
separately the average V2 for the triggers ((1*2)) and partners ((v%)). Due to eccen- 
tricity fluctuation, the vi of the triggers and partners fluctuate in the same direction 
event by event. This intrinsic fluctuation would leads to {v^v^) 7^ (w|)(wf). Both 
PHOBOSlSSland STARlSlfound a substantial eccentricity fluctuation, around 40%. 
From it, the additional correction factor can be roughly estimated to be (assuming 
a gauss shape for the fluctuation): 




1.08 (v 2 ) 



(4) 



So this suggest that the W2{RP} should be increased by 8% in order to account for 
the correlated fluctuation between triggers and partners. However, PHENIX have 
measured both the t>2{2} ^{RP}^^- The agreements between the two methods are 
better than 8%. Note that EqH]is very sensitive to the 5v2 used in the calculation 
and as well as the gauss smearing assumption. The correction factor would be only 
4.5% for a 30% input vi fluctuation. So probably the estimation of Eq0] is too 
simplistic. 

In many jet correlation analyses in STAR, the used vi is an average between 
TPC U2{RP} and «2{4}. Since the t>2{4} reduces non-flow as well as fluctuation 
effects, it is smaller than i>2{2}. This probably explains why their away side jet 
shape at intermediate pr is less concave comparing to PHENIX. 
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Fig. 13. The fake V2 of the leading particle from jet as function of centrality using the EP deter- 
mined in four rj windows. The embedded jet is at mid-rapidity (\rj\ <0.35). 
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One of the advantages of the PHENIX RP V2 measurement comes from the large 
pseudo-rapidity gap between BBC and the central arm. This gap greatly reduce the 
effects of the jet bias on the RP V2 ■ PHENIX have perform detailed Monte Carlo 
simulation to estimate the bias effect. Fig|T3]shows the centrality dependence of jet 
induced fake t>2 for various pseudo-rapidity window used to determine the RP. The 
fake V2 decreases as the r\ window used to determine the RP angle is further away 
from mid-rapidity. For the r\ window in BBC accept anc e (3 < |r/| < 4), the fake V2 
is almost negligible. Further details can be found in ES] 

PHENIX have measured the jet yield as function of angle with respect to the 
reaction plane at intermediate pr Due to the smearing effect due to relatively 
poor RP resolution of BBC and possibly a small jet signal, we see little jet like 
difference of the per-trigger yield in the in-plane direction relative to the out-plane 
direction(see Fig. 6 of ^ and discussions.). We can use this special property to 
constrain the V2 simultaneously using CFs measured at different angular direction. 
FigfT4k shows the CFs in 0-5% central Au+Au collisions for 6 angular bins in 15° 
steps. FigfMb shows the jet yield after subtracting the flow terms which can be 
calculated according to simple equation 1201. Although the CFs change dramati- 
cally from in plane to out of plane, the calculated flow term tracks the true flow 
background nicely. Given the small eccentricity in 0-5%, we can safely assume that 
the jet yield, jet yield should not depend on the trigger direction. Thus the rela- 
tively good agreements between 6 measurements is an independent confirmation of 
the t>2{RP} used in the correlation measurement. The small remaining difference 
between the jet functions in Fig[T4"b can be used to further constrain the V2- 
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Fig. 14. a) Correlation function for various 6 trigger direction bin and the trigger integrated bin 
(the center black curve) in 0-5% Au+Au at 200 GeV. b) The background subtracted per-trigger 
yields, the insert figure shows the 6 trigger bins. 



3.5. The meaning of per-trigger yield 

In two particle correlation analysis, typically one correlate particles in a high px 
window (type a) with those in a low pt window (type b). The distinction between 
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triggers and partners are arbitrary since medium modifications are on the jet pairs. 
There is a simple relation between the PTY using the high px particles as triggers 
and that using low px particles as triggers: 

RXaIIa = «aa4a = ^jt^et'patrspp = Ja ^Pt) (5) 

where JetPairsAA and JetPairs pp represent the average number of jet pairs in one 
A+A collision and one p+p collision, respectively. The i? AA represents Raa for 
type a particles and I AA represents PTY using type a particles as triggers. In 
naive jet quenching picture, the observed triggers are biased close to the surface, 
the corresponding away side companions have a larger medium to traverse. One 
expects /aa < -Raa or Jaa < RaaIaa- In direct photon-jet correlation, since 
direct photons are not modified by the medium, the away side jet would behave 
exactly as the single jet suppression: /aa = Raa- 

In correlation analysis, we are interested in studying the modification of the 
away side jet tagged by the triggering jet, i.e the per-jet yield. In reality we 
study the modification on the per-trigger yield, represented by Jaa- The physics 
interpretation of this quantity are complicated by the modifications on the triggers. 
In general there can be two types of modifications of the triggers: 

(1) Additional triggers not coming from jet: such as thermal-thermal recombina- 
tion, boost of soft particles due to flow. They tend to dilution the PTY. 

(2) Pure medium modifications of the triggers: jet quenching, enhancement due to 
energy dissipation such as radiative gluons/Mach-cone/Cherekov. They either 
does not affect or increase the PTY. 

(2) represents the effects we hope to study, which are complicated by (1). Both 
can change Iaa, but relation Eq[5]always holds. FigfT5lillustrates th idea by showing 
the PTYs from Au+Au and p+p, where the two hadrons are selected from non- 
overlapping pt ranges. In the left panel, 2-3 GeV/c hadrons and 3-4 GeV/c hadrons 
are used as triggers (a) and partners (b), respectively. In the right panel, the role 
of triggers and partners are swapped. Clearly, the vertical scale of the left panel is 
less than that in the right panel, simply because there are many more hadrons in 
2-3 GeV/c than 3-4 GeV/c. In the left panel, the nearside Au+Au data points are 
close to those for p+p; whereas in the right panel, the nearside Au+Au data are 
higher. This is the case because Raa is smaller in 3-4 GeV/c than in 2-3 GeV/c. 
More suppression means less number of triggers for fixed number of pairs, thus the 
PTY increase relative to the p+p value. 

Either (1) or (2) can explain this example. We can argue that most of type-a 
particles come from recombination without jet correlation, thus Jaa for type-a is 
smaller than that for type-b. We can also argue that type-a particles come from 
jet but are less suppressed than type-b, thus Jaa for type-b is larger than that for 
type-a. «2 and spectra measurements suggest that a large fraction (more than 50%) 
of the particles in 2-4 GeV/c come from recombination. However the measured 
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Fig. 15. The per-trigger yield for 2-3 and 3-4 GeV/c correlation, a) using 2-3 GeV/c particles as 
trigger, b) using 3-4 GeV/c particles as trigger. 



per-trigger yield is significantly enhanced at pt <4 GeV/c as shown by FigO This 
suggests that either triggers are dominated by the shower-thermal recombination, 
hence maintaining the jet correlation, or the enhancement of the jet multiplicity 
is much bigger than the dilution to the triggers. PID correlation measurement can 
help to clarify this issue. Fig[12p3 indeed shows some differences between baryon 
and meson triggers on the near side and possibly also on the away side. This trend 
is not consi sten t with the dilution from pure thermal-thermal recombination model 
predictions ^3. Further high statistics detailed measurement of the PID correlation 
in intermediate pr can help clarify the origin of the enhancement. 

4. Jet quenching 

4.1. hadron-hadron and tv° -hadron correlation at high pt 

At high pt, the situation is a lot more cleaner. The dilution to jet yield due to 
recombination or flow (if there is any) become negligible. Since each jet fragments 
into one high px trigger, the per-trigger-yield at the away side is the same as 
the per-jet yield. Fig[l6] show the centrality dependence of the jet yield for 5-10 
GeV/c trigger and 3-10 GeV/c partner in 200 GeV Au+Au collisions. The shape 
and yield at the near side show very little centrality dependence. The away side 
jet is clearly peaked at tt, but the magnitude is significantly suppressed in central 
Au+Au collisions comparing to that in peripheral collisions, consistent with strong 
quenching or absorption of the away side jet. 

PHENIX also carried out detailed study of the jet correlation analysis in Cu+Cu 
collisions at 200 GeV. The jet correlation measurements in Cu+Cu offer several 
unique advantages. In terms of iV part , Cu+Cu collisions cover from peripheral to 
30-40% central Au+Au collisions. Ap art can be determined with better precision 
that in Au+Au, thus allows a more detailed mapping of the centrality dependence 
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Fig. 16. The hadron - hadron PTY as function of centrality at high pt in 200 GeV Au+Au 
collisions. 



of the onset of the jet quenching. Secondly, systematic error from V2 on the jet yield 
is much smaller in Cu+Cu due to a smaller combinatoric background. 

The measurement was carried with 5-10 GeV/c ir° as triggers and charged 
hadrons in several ranges in 0.4-10 GeV/c as partners. FigfTTI show the correlation 
function in 0-20% Cu+Cu compared with that in p+p. A clear away side excess 
can be seen in all partner px ranges. No concave shape is expected even after the 
V2 background subtraction. At partner px > 2 GeV/c, the away side can be well 
fitted with a gauss function. Fig|18l shows the near side and away side jet width as 
function of partner pt for p+p, 0-20% and 20-40% Cu+Cu collisions. No apparent 
differences can be seen between Cu+Cu and p+p, suggesting partners come mostly 
from jet fragmentation in the px range under consideration (all px on the near side 
and px >2 GeV/c on the away side). 
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Fig. 17. The correlation function in Au+Au (upper panels) and p+p (bottom panels). 



Comparing the top and the bottom panels of Fig[T7j one notices that there 
is a larger asymmetry between the near and away side in Cu+Cu than in p+p. 
To quantify the medium modifications, we extract the PTY for Cu+Cu and p+p 
and construct the Iaa as function of xe {xe = Pt,a/Pt,t cos(A^i)). The results 
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Fig. 18. The near side and away side jet gauss width for p+p, 0-20% and 20-40% Cu+Cu collisions. 



from 0-20% central Cu+Cu collisions are shown in the left panel of FigOjJ] Iaa at 
the near side is consistent with 1 in the full Xe range of 0.1-1.4; the away side Iaa 
starts at slightly above 1 at small xe, and gradually decreases towards larger xe, at 
xe > 0.4 the suppression value approaches a constant of 0.5. This constant behavior 
is also seen by STAR in Au+Au collisions with A part ^4 The level of suppression 
was found to be similar to single particle Raa- Since the single particle Raa in 
Au+Au and Cu+Cu collisions follows iV part scaling, we can ask whether Iaa also 
has similar scaling behavior. The right panel of FigQj)] shows the integrated Iaa in 
0.4 < xe < 1 as function N part , comparing with the integrated Raa for the single 
particle. Assuming (z) ~ 0.7 for leading 7r°s, the original jets should be around 
5/0.7=7 GeV/c. Thus the Iaa of the away side should be directly comparable to 
the single particle at pr > 7 GeV/c. In reality, since the ir° Raa is flat at pr > 4 
GeV/c, the exact values of the jet energy for high px pi°s are not important. Near 
side Iaa is around one in all centralities, consistent with surface emission picture. 
Away side Iaa shows a suppression that has a similar centrality dependence in 
IVpart relative to single particle suppression. This is rather surprising given that 
the away side jet travels more medium than the single particles in the naive jet 
absorption picture. Probably suggesting that the simple geometrical bias argument 
in Sectior J3T5l is too naive. In jet quenching picture, the suppression is due to the 
energy degradation of the high px jets. The suppression factor -Raa depends on both 
the energy loss itself as well as on the input parton spectra shape. For the single 
particle, the expected iVbinary scaled p+p spectra have a typical power-law shape 
with a power of 8. In the di-hadron correlation, the away side spectra associated 
with the leading particles have much flatter distribution with a much smaller power. 
Thus for the same amount of energy loss for single jet and away side jet conditional 
to the near side jet, the suppression level observed for Iaa could be well less than 
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the Raa- If we follow the prescription in Ref. p , the fractional energy loss Si oss is 
related to Raa as: 

S loss = l-R A A(PT) 1/{n - 1} or Raa = (1 - S^Y 1 - 1 (6) 

The power "n" is 7.1 for single particle spectra in dN/dpr- Th e power for the away 
side associated spectra can be determined from p + p data^and it is about 4.8 in 
dN/dpr - If the away side jet have same fraction energy loss as single particle, then 
we would expect: 



'los; 



1 - Raa (pt 



,i/(«k-i) 



1 - Iaa{pt 



(nr-l)/(n K -l) 
AA 



(7) 
= 0.37, 



-Raa = 0.2 in central Au+Au collisions would lead to Iaa = R 
much bigger than Raa- On the other hand, if we require Iaa — Raa, then the away 
hadron energy loss fraction would be S{ oss = 1 — /aa(pt) 1 ^™ 1-1 ^ = 0.345, much 
bigger than the single hadron energy loss fraction S^ BS — 0.23, as expected (about 
50% more energy loss). If we apply the same trick to central Cu+Cu collisions, 
assuming Iaa — Raa — 0.5, the fractional energy loss is would be S^ ss = 0.107 
and S( oss — 0.167 for single particle and away side jet, respectively. 



4.2. j-hadron correlation 

In high pt hadron - hadron or 7r°-hadron correlation, the total jet energy is not 
known. Trigger hadron from normal jet fragmentation carries on average about 60- 
70% ((z)) of the total jet energy. This effect is known as the trigger bias, which 
have been studied in detail in p+p and d+Au collisions 1 ^9 1 17 | ener gy \ osa f 
triggering jet in the medium leads to additional bias effects, and since the energy 
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loss is dependent on path length, the trigger bias effects would also coupled with the 
collision geometry. Thus complicates the extraction of the medium effects. Direct 
photon (7dir)- hadron correlation does not suffer from such limitations. 7di r serves 
as the gauge of the jet energy and jet direction, allowing the direct study of the 
medium modification of the away side jets. In reality, 7di r — h correlation has its 
own difficulties and limitations. In particular, once have to deal with large decay 
7 background, non-hard-scattering contributions such as the fragmentation and 
bremsstrahlung and thermal radiation and small rate. 

PHENIX has observed a pt dependent direct photon excess above the decay 
background in central Au+Au collisions. A sta tistical subtraction method combined 
with the knowledge on the direct 7 excess is used to obtain the 7^ — h signal. 
Figi20]show the hadron yield associated with direct photons in p+p collisions for two 
partner pt ranges. Near side is consistent with zero, and away side shows some finite 
excess. The p+p data are compared with the expected values from the PYTHIA6.1 
calculation with fcy=2.5 GeV/c. Within the large systematic uncertainties, the data 
qualitatively agrees with the PYTHIA value. 
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Fig. 20. The direct 7-h signal in p+p collisions compared with PYTHIA simulation for 7-9 GeV/c 
trigger 7 and two partner pt ranges. 



Left panel of FigJ5Tp21 shows the comparison of the A<f> distribution of the PTY 
in p+p and 0-20% central Au+Au collisions. The Au+Au data have smaller system- 
atic errors, however the measurements are currently statistics limited, especially at 
around A0 = ir/2 region. The Au+Au data seems to show a small positive yield at 
7r, but the significance is only about 1 a level. One can reduce the statistic errors 
by integrating the away side yield in ir/2 < Acj) < 3ir/2. The right panel of Fig|2"T1 
show the integrated yield for several trigger px regions with partner px fixed in 2-5 
GeV/c. p+p data suggest a gradual increase of the PTY towards larger trigger pt, 
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Fig. 21. a) The per-trigger yield for direct 7 in central Au+Au (circles) and p+p (boxes), b) The 
integrated yield of the away side for several trigger direct 7 ranges (data are plotted in the middle 
of the bin) 



while the Au+Au data seems to stay relative flat and below the p+p level. 

Overall, the results of the direct 7 — h correlation is still in a very qualitative 
stage, mainly due to limited statistics. However, the necessary methodologies have 
been developed and tested by PHENIX. It is only a question whether there will be 
sufficient statistics from future Au+Au runs. We might need the RHIC-II in order 
to do a really jet tomography measurement with direct photons. 



5. Conclusion and outlook 

PHENIX have mapped out the landscape of two particle azimuth correlation in 
heavy ion collisions as function of pt,t, Pt,a, centrality, hadron species, \fs and 
collision systems in different A</> regions. The background subtracted "jet" pairs 
seems to come from four different components concentrated in different Ac/> regions 
with their characteristic dependence on px, PID, y/s etc. 1) A hard component 
around A(j> = that is consistent with jet fragmentation; 2) A soft and broad 
component around A</> = that is consistent with the "Ridge" seen by STAR; 
3) A hard component ("Head") around Ac/) = tt consistent with fragmentation of 
away side jet; and 4) a soft component centered at \A<j> — tt| ~ 1 ("Shoulder"). The 
hard components are sensitive to the "quenching" of the (di-)jets by the medium. 
The hard component at the near side shows little modification, whereas the hard 
component at the away side are strongly suppressed, up to factor of 5 in central 
Au+Au collisions. The soft components ("Ridge" and "Shoulder") represent the 
response of the medium to the jets, appearing as distortions of the shape and 
enhancements of the yield. These components are shown to be important at px < 4 
GeV/c and have a particle composition that is closer to the bulk medium. The 
non-trivial evolutions of the near and away side shape/yield with px and y/s are 
the results of the detailed interplay between the soft and hard components. 

The four different components can have very different geometrical bias and 
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trigger bias. In a simple jet absorption picture where a very opaque medium is 
assumed, the near side hard component is emitted from the surface, (L) ~ 0; the 
near side "Ridge" comes from jet close to the surface, (L) < R (R is the average 
radius of the medium); the away side "Shoulder" results from jet that traverses a 
large path length, (L) > R; the path length for away side hard component can vary 
dramatically depends whether it is tangential emission or punch-through, < (L) < 
2R. One of the important future directions of the jet correlation is to quantitatively 
separate the four different components; and to study the detailed properties of each 
components. To fully understand the medium response, one need to address the 
question that whether the near side "Ridge" and away side "Shoulder" are of the 
same origin or not. A detailed mapping of the px, PID, charge and y/s dependence 
would be very helpful in this regard. One can dial the path length by triggering on 
two high pr triggers and correlate with the third soft hadrons, the "displaced" peak 
might be seen on both the near and away side SD. However this is only possible if 
the high pt have a significant punch-through component and the medium can not 
be very opaque. On the jet quenching part, high pt correlation remains to be a good 
tomographic tool since it is less affected by the surface bias Both the centrality 
and reaction plane dependence of the jet correlation at high px would be very 
helpful in constraining the transport properties of the medium. In the meanwhile, 
we should continue to pursue the gamma-jet correlation with increased statistics 
and refined analysis techniques. Since it is less affected by the geometrical bias and 
trigger bias, even a statistics lim ited result could be very powerful in constraining 
various jet quenching models EH. 
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